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Abstract 

We describe the "Feynman diagram" approach to nonrelativistic quantum mechanics on R" , 
with magnetic and potential terms. In particular, for each classical path 7 connecting points 
go and qi in time t, we define a formal power series Vy{t,qo,qi) in fi., given combinatorially by 
a sum of diagrams that each represent finite-dimensional convergent integrals. We prove that 
exp(l/y) satisfies Schrodinger's equation, and explain in what sense the i — > limit approaches 
the i5 distribution. As such, our construction gives explicitly the full fi. — > asymptotics of 
the fundamental solution to Schrodinger's equation in terms of solutions to the corresponding 
classical system. These results justify the heuristic expansion of Fcynman's path integral in 
diagrams. 
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1 Introduction 

The goal of this paper is to show that the semiclassical description of the evolution operator in 
quantum mechanics in terms of Feynman diagrams satisfies all the necessary requirements: it sat- 
isfies the Schodinger equation and has the correct initial conditions. This description is rooted in 
Feynman's path integral interpretation of quantum mechanics, first published in [■'5] (see also [G]), 
and formal semiclassical arguments. In certain cases, non-perturbative path integrals can be inter- 
preted analytically as Weiner integrals, e.g. [1, 7, 10]. However, for the last half century, physicists 
have used "formal" perturbative path integrals and Feynman diagrams extensively to make per- 
turbative empirical predictions in theories that do not have adequate mathematical foundations, 
and there have been similar mathematical applications to the subject of topological quantum field 
theory [15]. Surprisingly, the diagrammatic methods have not previously been shown to give the 
correct answers when applied to the one mathematically rigorous nonperturbative theory available, 
namely nonrelativistic quantum mechanics. We fill that gap. 

In this paper we consider quantum mechanics in M" with the Hamiltonian being a second-order 
differential operator. (In a [>-] we discuss perturbative quantum mechanics on general Riemannian 
manifolds.) Recall some basic definitions. Let Bi, . . . , Bn,C : — )• M be smooth functions.^ 
The nonrelativistic Schrodinger operator with electric potential C and magnetic potential B = 



^Without much difficulty, these functions can be taken to depend explicitly on an external time parameter, but 
for notational convenience we do not do so. The function U{t,qo,qi) that we will discuss should then be a function 
of two time variables U{to,ti,qo,qi), and similar superficial changes must be made throughout this paper. 
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{Bi, . . . ,Bn) is the following linear second-order differential operator on M": 

Solutions to the partial differential equation ^ip{t,q) = jf^H'il'{t,q) with fixed initial data define a 
one-parameter family of unitary "time evolution" operators U{t) on L^(M"). The kernel of the time 
evolution operator is a distribution U{t,qQ,qi), satisfying {U{t)ijj){qi) = Jjg„ U{t,qQ,qi)ilj{qQ)dqQ. It 
is determined by the following initial value problem: 

ifi^^U{t,qo,qi) = Hg,[U{t,qo,qi)] U{0,qQ,q) = 5{q - qo) (1.0.2) 

The operator Hg-^ acts on the gi variable, leaving (t, qo) fixed. In this paper we describe explicitly 
the asymptotics or U{t, qo, qi) as ^ — )• 0, at least for (t, qo, qi) in a dense open subset of M>o x M^"'. 

The asymptotics of U{t,qo,qi) are known to be given by the Hamilton-Jacobi action and solu- 
tions to the transport equation [16]. Feynman's path integral also predicts that the asymptotics 
are given by a power series where coefficients are parameterized by "Feynman diagrams": each 
diagram represents a finite-dimensional integral, and the coefficients of the power series are the 
values of these intergals. In this paper we will define this diagrammatic expansion; the integrands 
depend on the solutions to the corresponding classical mechanical system. We will prove that the 
Feynman-diagrammatic power series satisfies formal versions of equation (1.0.2): it satisfies the 
Schrodinger equation as a function valued in formal power series, and as a distribution valued in 
formal power series approaches the 5 distribution as t — )• 0. 

In select cases these results are already known. For the quadratic "harmonic oscillator" (the 
functions Bi,. . . , Bn are linear in position, C is quadratic), the path integral is defined directly from 
Gauss's formula, and agrees with the Schrodinger picture [ ] . Feynman's original argument applies 
when the corresponding classical physics is a perturbation of the free theory, i.e. when Bi, . . . , Bn, C 
are infinitesimal quantities. Combining these observations, it is not too difficult to prove that the 
diagrammatic path integral satisfies Schrodinger's equation for infinitesimal perturbations of the 
harmonic oscillator. More importantly, Duru and Kleinert [ :] apply the path integral approach to 
calculate the energy spectrum of a particle moving in the Coulomb potential, and compare their 
approach to the Schrodinger methods. Their methods differ from ours: we disallow singular po- 
tentials but consider arbitrary smooth ones, so our results do not apply to the Coulomb potential; 
moreover, they take advantage of numerous changes of variables (without explicitly checking the 
coordinate invariance of the path integral) and use the "charge" of the potential well as the per- 
turbation parameter, whereas we define the path integral directly and use Planck's constant h for 
perturbation. Diagrammatic path integrals do exist in the work of Kleinert and Chervyakov [1 1] 
and in the work of DeWitt-Morette [-]; Kleinert and Chervyakov's methods apply directly only to 
perturbations of free motion, and neither they nor DeWitt-Morette (who does give an expansion 
equivalent to ours) check the Schrodinger equation directly, although the result is implicit in their 
and Duru and Kleinert 's approaches. 

We summarize the main definitions and results in Section 2. Section 3 has proofs of two key 
lemmas. The main theorem is proven in Section 4, and in Section 5 we prove the short-time 
asymptotics. 
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2 Statements of the definitions and results 
2.1 Preliminaries 

We begin by establishing some notation. Kronecker's 5-matrix is 6j = liii = j and otherwise. 
Dirac's (5-function is the distribution defined by f {x)5{x)dx = /(O). Heaviside's step function 



-oo 

is: 



e(x) = / 5{t) dt 



0, X < 

1, x = 
1, x>0 

We henceforth adopt Einstein's summation convention: piq^ = q^pi = Yll^=iPiQ^- We raise and 
lower indices as is convenient: = Pi- 

2.1.1 Definition A (piecewise-smooth parameterized) path is a continuous function 7 : 
[0,t] —7- M" such that there exists a decomposition = to < ti < ■ ■ ■ < tn = t with jlit^^tj+i] smooth 
for each j = 0, . . . , n — 1. The duration of a path 7 : [0, t] — )■ M" is the real number t. A based 
loop is a path 7 : [0, t] — )• M" with 7(0) = = "f{t). A path 7 is classical if it satisfies the following 
nonlinear second-order differential equation (we write 7(t) = ^7(t), etc.): 

Let 7, ^ he paths of duration t, and define the differential operator given by: 

A classical path 7 of duration t is nondegenerate if the kernel ofDj among based loops of duration 
t is trivial. 

As equation (2.1.2) is nondegenerate second-order, a classical path 7 is determined by its du- 
ration and its initial conditions (7(0), 7(0)). Thus the space of classical paths is naturally an open 
subset of M2"+i_ Nonde generacy is clearly an open condition among classical paths. A standard 
result holds that for given initial conditions, the set of durations for which a classical path fails to 
be nondegenerate is discrete; see e.g. [I ]. In particular, the nondegenerate classical paths are a 
dense open subset of all classical paths. The following lemma follows from e.g. [ )]; in Section 3.1 
we include the proof rewritten in the language of classical paths. 

2.1.4 Lemma Letj : [0,t] — )• M" be classical and nondegenerate. Then there exists a neighborhood 
O of (7(0), 7(t)) G M^" and a map 7 : [0, t] x O — ;> M" such that (i) 'j{—,qo,qi) is classical 
and nondegenerate for each {qo,qi) G O, (ii) j{0,qo,qi) = qo and ^{t,qo,qi) = qi, and (Hi) 
7(-,7(0),7(t)) =7. 
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Let 7 be a classical path of duration t. Then for sufficiently small e, 7 determines classical 
paths 7<j for s G [t — €,t + e], either by restriction or by extension via equation (2.1.2), and if 7 
is nondegenerate, for sufficiently small e so are the 75. Thus, to a nondegenerate classical path 7 
we can associate a family itiQojQi,—) of nondegenerate classical paths, varying smoothly in the 
parameters. We will abuse notation and write 7 for this whole family. So henceforth 7 is a smooth 
family of solutions to the boundary value problem given by equation (2.1.2) along with Dirichlet 
boundary values 7(0) = qq, ^{t) = qi, where {t, qo, qi) ranges over an open subset O of R2"+-'^. 

2.1.5 Definition Let j be a classical nondegenerate path, extended to a family with boundary 
values varying in O CI M?"'~^^ as above. The corresponding Hamilton function S*^ : O — t- M is: 

S^(.t,qo,qi) = 1^ Q|7(r)|' +i?.(7(T)) Y{t) - C(7(r))) dr (2.1.6) 
Then Sy{t,qQ,qi) satisfies the following well-known equations: 

--l\Vg,[-S,]-B{qo)f + C{qo) = l\V,d-S,]+B{q,)f + C{qi) (2.1.7) 



dt 



dq'o dq\ 

For a few different constructions, we will need the matrix q^q^^ [—S{t, qQ, gi)] and its inverse; we 

\ dqidqa ) J dq' dq"^ ~ similarly for qi ^ qo- 

2.1.9 Lemma Let j be a nondegenerate path of duration t, extended to a family as above. Let 
: [0,t]^2 ^ (i^n)®2 defined by: 



dl\ ^/^/52(-S'^)^-lY'a7 



dqo V- dqodqi J J dq[ 

Then —G^ is a Green's function for the operator T>^ given in equation (2.1.3): i.e. G^ satisfies 
P-y[G^]^(?, r) = —5\.5{t — q), where T)^ acts on thej,T parts of equation (2.1.10), and G^ vanishes 
along the boundary of the square [0,t]^'^. 

For the proof, see Section 3.2. Statements similar to lemmas 2.1.4 and 2.1.9 appear in e.g. [2], 
where degenerate paths are also briefly addressed. 



2.2 Feynman diagrams and the main results 

2.2.1 Definition A partition of a set X is a collection of nonempty subsets of X, called the 
blocks of the partition, that are pairwise disjoint and whose union is all of X. 

A Feynman diagram is a finite combinatorial graph with all vertices trivalent and higher. 
(Self-loops, parallel edges, and the empty diagram are allowed.) More precisely, a Feynman diagram 
is a finite set H of "half-edges", along with two partitions E (the "edges") and V (the "vertices") 
of H , with the requirements \e\ = 2 for each block e £ E, and \v\ > 3 for each block v €z V . 
An isomorphism of Feynman diagrams Ti = (Hi, Ei,Vi) and T2 = (-^2, -E'2, V2) is a bijection 
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(j) : Hi ^ H2 such that 4>{Ei) = E2 and (l){Vi) = V2, where (j){Ei) , cj){Vi) are the obvious partitions 
of H2 induced by {(j),Ei,Vi). We will denote Feynman diagrams with the obvious pictures: 



Let r = [H, E, V) be a Feynman diagram. The Euler characteristic ofT is x(r) = |^| — \E\. 
The connected components of T are the equivalence classes of H under the equivalence relation 
generated by EUV . Let ttq^T) be the set of connected components ofV; then A(r) = |7ro(r)| — x(r) 
is the number of loops of T , also called its first Betti number. A diagram T is connected if 
it has precisely one connected component. For example, the above three pictures comprise all (up 
to isomorphism) diagrams T with = —1, and each is connected. 

A marked Feynman diagram is a Feynman diagram T = (H, E, V) along with a subset 
M Q H of "marked half-edges" such at each vertex v £ V , at most one half edge rj £ v is marked, 
i.e. \vriM\ < 1. An isomorphism of marked Feynman diagrams is required to respect the marking. 
We indicate the markings on a marked Feynman diagram with small circles, as in equation (2.2.4). 

2.2.2 Definition Pick a duration t > 0. A totally labeled Feynman diagram (with duration 
t) is a Feynman diagram T = {H, E, V) along with maps t : V ^ [0, t] and j:H^{l,...,n}. Let 
^ be a classical nondegenerate path of duration t and (r,M, r, J) a totally labeled marked Feynman 
diagram. Then its value is (no summation convention): 



J-^(r,M,r,j)= n 

v€V, i;nM=0 



X n 

vGV,vnM={(} 




+ n 



d 



[Bk{q)] 



j^T) I X 



g='r{r{v)) 



\ \vev~.{<:} 



enM=0 



e={r,,C}eE, 
enM={ri} 



n 



?2-t(C) gCM 



dc^idq2 



?i=t'(»?) 
f2=T{C) 



//(r, M) is an (unlabeled) marked Feynman diagram, its value is an integral over all total labelings: 

dr 



Je{l,...,n}^^ ^g[o,t]V 



More compactly, Feynman diagrams are evaluated via the following Feynman rules and the ob- 
vious compositions and contractions (including the Einstein summation convention), where 5,i,...,ik 
are test-functions ranging over the space of paths of duration t: 



d^Ciq) ^ d^Bkjq) 



\.^Q \dq^^ ■ ■ ■ dq^" dq^^ ■ ■ ■ dq^" 

66 X o„-l 



i\'{r) ■■■Cl,Hr)dT (2.2.3) 



d^-^[-Bi,] 



T=o dq^^ ■ ■ ■ dq^" 



9=7(t) 



(2.2.4) 
(2.2.5) 
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We must check that the value of a marked Feynman diagram converges: the Green's function 
r) is non-smooth Hke |? — t| near the diagonal, so any diagram in which both ends of the same 
edge are marked includes a Dirac 5- function distribution. But no vertex can have more than one 
marking, so such (5-functions do not lead to divergence integrals. (We have insisted that B, C be 
smooth. If they are allowed to be singular, divergences can occur, requiring some renormalization 
procedure, c.f. [12].) We prove the following in Section 4. 

2.2.6 Theorem (Schrodinger's Equation) Let j be a classical nondegenerate path, extended to 
a family with boundary values {t,qQ,qi) ranging over some open set O G M^"+-'^. Define on O the 
following M\[ifi\\- valued function: 



V^{t, qo, qi) = -S^{t, qo, qi) + ih^ log 



det 



d'\-s. 
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dqodqi 



E 



, , |Aut(r,M)| 

equivalence classes (r,M) 
of connected marked 
Feynman diagrams 

(2.2.7) 

Then U-y = exp[{ifi)~^V^^ is a formal solution to Schrodinger's equation: ih^^ = H[U^] as formal 
power series in h. 

We remark that for a connected diagram T = {H, E, V), we have A(r) = \E\ — \V\ + 1 < \E\ = 
2\H\, and so the sum in equation (2.2.7) is finite at each order in h. In general, the series has zero 
radius of convergence in h. Note that U-y can be described as: 



Uy{t,qo,qi) = e 



-(ih)- 




det 



d[-Sy] 



dqodqi 



equivalence classes of 
marked Feynman 
diagrams {T,M) 



|Aut(r,M)| 



2.3 Sums and limits of Uy 



Temporarily, let h range over "small" non-zero real numbers. Pick a positive integer and suppose 
that O C M" is open, and /i,/2 : O ^ M.[ih] satisfy /i = + 0{h^+^). Let : O ^ R" have 
compact support. Then: 

j f,m-'h{q) g^q) dq = j e^'^^''^'^"^ g{q) dq + 0(?i^+"/2) (2.3.1) 

For details, see [4, Chapter 3]. Thus, if f,g : O ^ M.[[ih]] are smooth term-by-term and g has 
compact support, then we are justified in defining f ex.p[{ih)~^ f{Q)) g{<l) dq as a formal expression 
in h, by successively approximating / by polynomials in h. We will write J*^"^™^^ for this formal 
integral. 

We will not totally define what we mean by a "formal expression in although it is not hard 
to do so. Essentially, a formal expression in h is an C[[-v/^]] -linear combination of terms of the form 
exp((i?i)~-'^/) , where / G M[[i/i]], modulo the obvious equivalences and with the obvious arithmetic. 
If / (mod h) has finitely many critical points within the support of g, and if at each critical point 
the Hessian of / is nondegenerate, then exp((i/i)~^/) g is a. formal expression in h in this 

sense. 



2.3.2 Definition Let /i, /2, . . . , foo be distributions on O Q R", valued in formal expressions in h 
Then lim /jv = foo if and only if for every g : O ^ M. with compact support, limAr^oo J^™™'^^ fNiq) gig) dq 
J^'^^^^ foo{x) g{x) dx. This defines the pointwise convergence topology on the space of distribu- 
tions. The topology on W[[ih]] is given coefficient-by-coefficient. 
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In Section 5 we will prove: 

2.3.3 Theorem (The initial value problem) Pick O C r2"+^ open, and let A he the set of 

all families of nondegenerate classical paths with boundary values varying over O. Let {g-y}-ieA 
be a collection of smooth functions O — >• ]R[[i/i]]. Then for fixed t,qi, the possibly-infinite sum 
^^j^Adiitjqojqi) exp((i?i)^^ V^(t, go, gi)) converges pointwise as a distribution in qo valued in for- 
mal expressions in h. 

Moreover, suppose that O = (0, e) x Oq x Oq, where Oq C R" is open with compact closure 
and and e > is sufficiently small depending on Oq . Then for fixed qi, as a distribution valued in 
formal expressions in h we have: 

lim V exp{iih)-%it, qo, qi)) = {27ihT/\'^^'%qo - qi) (2.3.4) 

The reader should be mildly disappointed in Theorem 2.3.3: we have first set ^ to a formal 
parameter, and then, in equation (2.3.4), taken a limit as t — )• 0. Much better would be to approx- 
imate by a polynomial in h, interpret ^ as a positive real parameter, take any necessary limits, 
and then take the asymptotics as — >• 0. But when the limits are performed in this order, we 
do not believe that X^^g^i exp((i/i)~^y^) always converges. Moreover, even if it does, some extra 
assumption is needed in order to take the i — )• limit for non-formal h in equation (2.3.4). For 
example, it is enough if the paths 7 G A are of bounded length, although we will not prove this. 
However, in as natural a system as the one-dimensional mechanics with electric potential C{q) = q^, 
the classical paths can be arbitrarily long, and as t — )• with fixed endpoints the classical paths 
run away with length growing as t"-*^/^. 

One final remark is in order. It is tempting to conclude from Theorems 2.2.6 and 2.3.3 that the 
^ — )• asymptotics of the full propagator U (t, qo, qi) are be given by: 

U{t,qo,qi) « i27Thr^/h^^^/^exp{iihr'V^,{t,qo,qi)) 

7 classical, 
7(0)=go,7{*)='?i 

But a modification must be made. The formal power series are defined only for 7 nondegenerate, 
and Theorem 2.2.6 only applies in this case. Thus the above formal sum solves equation (1.0.2) 
only for short times in a neighborhood of the diagonal. In general, if we know the asymptotics of U 
near short nondegenerate paths, we can determine the full asymptotics by means of the composition 
law: U {to + ti,qo, qi) = J U{to, qoiQ) U{ti, q, qi) dq. In [n] we prove the following: 

2.3.5 Theorem (Semigroup Law) Let 7,V^ be as in Theorem 2.2.6. Suppose that t = to + 
ti and the restrictions 70 = 7|[o,to] '^'^^ 7i — l\[toA 0.''^^ nondegenerate, and define the corre- 
sponding formal power series and V^-^ . Let 77(7) be the Morse index of 7 with respect to 
the Morse function 7 1— t- /q(^|7|^ + -6(7) • 7 — C{j)^dT on the space of paths with the same 
boundary values, and similarly define Morse indexes r]{'yo) o-^d r/(7i) (see e.g. [ ]). Set [/^ = 
(_^)f?(7)(27r/i)-"/2e~™'^/*^exp((i^)~-^ y^(t,qo,^i)); ^.n-d similarly define U^y^ and U^^. Then for a 
sufficiently small neighborhood Q 0/7(^9) we have: 

/■formal 

U-y{to + ti,qo,qi) = / U^oito,qo,q)Uyiiti,q,qi)dq 
JQ 

Thus, the correct asymptotics of the full propagator U{t,qo,qi) are essentially given by: 
Uit,qo,qi) « (27rn)-"/2e^"-/4 ^ (_i)r,(7) exp((in)-i V,{t,qo,qi)) 

7 classical, 
7(0)=iJ0,7W=9i 
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Here the sum makes sense provided that there is no degenerate classical path of duration t con- 
necting qq to qi, and converges pointwise as a distribution. This justifies the method of formal 
Feynman path integration. 



3 Proofs of the Lemmas 
3.1 Proof of Lemma 2.1.4 

Let V be the space of paths of duration t and A the space of based loops of duration t. Each is an 
infinite-dimensional vector space; we think of V as an infinite-dimensional smooth manifold with 
tangent bundle TV = V x V. Recall that the derivative df of a function / : P — t- M is defined 
by d/7 • C = ^1^=0 ~'~ ^ differentiable if such a derivative exists for all (^,7) G TV. Let 
TT : "P — )• M^" be the map 7 1— )• (7(0), 7(t)); it is a trivial bundle with fiber A. As A is a vector space, 
we can identify the fiber tangent spaces: if 7 G then T^(7r-^(7r(7))) = A. 
Consider the function / : P — )• R given by: 

fil) = 1^ (^|7(t)|' + i?.(7(T)) f (r) - C(7(r))) dr 

As is well-known, a path 7 is classical if and only if A C kerd/^. The Hessian of / at 7 is the 
symmetric bilinear form: 

: (6,6) ^ 1^ (^i(T)^^(r) + {i{{r) Ciir) + e^r) Ch{r)) + 

When 6)6 ^-re based loops, integration by parts gives h~f{^i,^2) = Jq^i{'t) " '^'y[^2]{T) dr . Thus 7 
is nondegenerate if and only if A n ker h-y = 0. 

Let C be the space of classical paths. By considering the derivatives of df along paths in C, one 
easily checks that T^C C ker h^. Thus if 7 is classical and nondegenerate, then dir : T^C — )• T^j-^^R^" 
is an injection. But the space of classical paths of a given duration has dimension 2n, and so dir is 
a bijection. This completes the proof of lemma 2.1.4. 



3.2 Proof of Lemma 2.1.9 

We begin with the observation that for each o = 0, 1 and each j = {1, . . . ,n}, the path (pajiT) = 
^j{t) satisfies P'y[i?5>] = 0; this follows by differentiating equation (2.1.2). In other words, the 

2n paths (paj are Jacobi fields along 7. The boundary values are: (pQjiO) = 6j = (j)\j{t) and 

^hj(^) ~ ~ Let 4>o{t), (t)i{T) be the n x n matrices {(f'ai'T'))^- = <P\j{T)- Each solution 

to the initial value problem 9^(0) = 0} is determined by 0(0). Since 0i(t) = 5 is full-rank, 

01 (0) must be also. 

Define the following function on [0, t] valued in 2n x 2n matrices: 



'^^"^=^-(r) Mr) 



Then M(0) is lower-triangular with full-rank blocks on the diagonal, and hence invertible. Since 
the coefficient on ^ in P^f^] is traceless, Liouville's formula gives det M(r) = detM(O), and so 
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M(r) is invertible for each r E [0,t]. Let ipo,ipi be the (n x n)-matrix valued functions on [0,t] 
satisfying: 



Then it is straightforward to check that the Green's function for —D^ is given by: 
Gii,, r) = -e(? - r) <pyr) + Q{t - ,) 0]'^(r) ^^1,(0 

But the differential operator —V^ is self-adjoint for the inner product (^1,^2) = /q " on 
the space of based loops. Thus if G\{q,T) is its Green's function, then GJ,(<j,r) = G^(r, In 
particular, each if^a satisfies ©^[V'a] = 0, and so the V'aS are linear combinations of the ^^s. To 
find the coefficients, we need only check the boundary values. We have ipQi^iO) = = V'j ^(i), by 
evaluating the top row of equation (3.2.1) at r = 0,t. On the other hand, the bottom row then 
gives: 



{Mt))-' = Mt) - ^'^^^^ 



drdqo 



T=t 



.-(m = ^ - = (3.2.2) 

dqo dqo \ dqi J dqodqi 



We have used equation (2.1.8) to do the necessary simplifications. Thus -00 ki'^) — s^i jv'^ ; y\ dqodq 

A similar argument shows that '4^\^{s) = ~'Poi{'^) ii dq^dqo ^ j ' "^^^^ completes the proof of 
lemma 2.1.9. 
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4 Proof of the Main Theorem 

The Schrodinger equation for exp((i/i)~^Ky) is equivalent to the following nonlinear differential 
equation for Vy-. 

dV.y 1,„ .^^,0 ih . . dVy ihdBj(qi) /1,„, ^.o A 

Our strategy will be to convert equation (4.0.1) into the diagrammatic language introduced earlier, 
whence it will follow from a few calculations of derivatives and some simple combinatorics. 

4.1 More Feynman rules 

In addition to the Feynman rules given in equations (2.2.4) to (2.2.5), we introduce the following 
notation to our graphical calculus: 

i 

) = 6'^ -^^((g)) = C{qi) T,(^(b)^ = BM 7 ) = l\r) 

i j I 

T,i 

In particular, dashed lines carry only an index (no time variable). Contraction of indices is implied 
by connecting dashed strands. 

We will henceforth largely drop the "J-".y" notation from equations of diagrams. We will find it 
convenient to abbreviate the two kinds of vertices in marked Feynman diagrams: 



9 



Then rather than discussing marked Feynman diagrams, we will simply say "Feynman diagram" 
and mean a diagram drawn with a-type vertices. It's clear that the automorphism counts work 
out: if adding some markings to a diagram divides the size of the autmorphism group by some 
number n, then there were n equivalent ways to add those markings to the unmarked diagram. 
In particular a symmetry of a diagram with a B-type vertex either acts as a symmetry on each of 
the expanded diagrams or permutes the possible expansions, and those that it permutes have their 
symmetry groups broken by the right amount. 

We will also extend the possible valence of such vertices in our notation, although not in the 
sums of diagrams: 



[ = (f (0) + BM)e{o) - {fit) + BM)e{t) 

+ 



+ 



e(T)-C(r)dr 

Here the test functions ^, ^ are paths of duration t, and in the second line on the right-hand side 
the bivalent vertices have the natural meanings extending equations (2.2.3) and (2.2.4) to n = 2. 

In particular, ^ = if ^ is a based loop. By integrating by parts, \^ = /q C(t) • P^[C](T)dr 
whenever ^ is a based loop, where T>^ is the differential operator given in equation (2.1.3). 
We represent differentiation with respect to t,qo,qi with dotted circles: 




d_ 

di 



Then the product rule can be written graphically as: 




d 



+ 



Suppose then that Fi is a subdiagram of F whose images Fi, . . . , F„ under the group of automor- 
phisms of F do not intersect, so that F = F U Fi U • • • U F„ (we do allow automorphism of F to 
induce nontrivial automorphisms of Fi). Then: 



; F ; = ■., F Fi F2 • • • Tn J = ■ F^Fi F2 • • • F„ + Ff Fi iFs • • • F„ + • • • + F Fi F2 

= (r) Fi F2 • • • F„ + n f (Fi) F2 • • • F„ 

It is an elementary counting lemma that Aut Fi ■F2 • • • F„^ = ^ Aut (FFi F2 • • • F„). From 
this observation, we derive the fundamental result that for z = q\ or t: 



- y 



dz ^ lAutFl 

r connected 



r connected with 



lAutFl 



(4.1.1) 



one .'iy-, or > 



Pi 
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The sum on the right-hand side ranges over connected components with precisely one differenti- 
ated basic subgraph — either a single differentiated vertex (of valence three or more) or a single 
differentiated edge. 

In this notation, equation (2.1.7) reads: 



1 Ql.' "-^H 91^' 
:T :■} 



(4.1.2) 



Moreover, it's clear that = ~^^7*(*)- Given equation (2.1.8), we have: 



^ ^ ' ~ - ^ (ji 91 ^ ' 



(4.1.3) 



Finally, we introduce the Feynman rule o'C^^JJ^\~ for the inverse matrix 



Q-T— — 



qi y.::l qo 



m-s])- 

07—- ri 



qi y---.f,qo 



For example, in this notation equation (2.1.10) reads: 



,.90 91 J. Q((; - r) + J. 91 



{d'[-S])-' 



*'.J. Q{t - ?) 



(4.1.4) 



4.2 Derivatives of the vertices 

In this section we describe the derivatives of the vertices. We begin with the derivatives with 
respect to qa- Equation (2.1.8) implies: 



9a 



:7 :■ 



Integration by parts and the observations in the first paragraph of Section 3.2 give: 



9a ^ 




dqa 



T=0 



In valence n > 2, the derivatives are particularly easy. Differentiating under the integral sign 
and using the chain rule gives: 



d d 



q='y{T) dqa 



.1 ?1 ?2 Cn 

dr= (Ti- 



ll 



A similar argument gives: 




Summing and noting also that the non-diagrammatic part of the bivalent vertex does not depend 
on the QaS, we have, for all n: 




(4.2.1) 



The derivatives with respect to t are harder because of the boundary terms. At valence n = 0, 
equation (4.1.2) already describes * f ■ }. We will not describe the other t-derivatives in full generality, 
as we will only need them for valence n > 3 when all incoming edges are based loops. In this case, 
when n > 3 at least one incoming edge is not marked, and hence the boundary terms vanish. We 
have: 

?1 f 2 t . ^1 ^2 in 

^S^^^ = ) n>'i, all ^jS are based loops (4.2.2) 

We conclude this section by addressing the second derivative of a vertex with respect to the qa- 
By the product rule: 






In particular: 




:7:- 7 : 



+ ■■: 



;7 




if ^ is a based loop. 



The Green's function G^, which we denote by an edge, is a based loop in each variable, so: 



.1. J. 

We have observed already that when contracted with based loops, the bivalent vertex acts as P-y. 
But 2?^[G-y]* (?, r) = — (5*(5(<j — r) by definition. In short-hand: 

« { 

, provided ^ is a based loop. 



(4.2.3) 



And^ 



vanishes at both endpoints r = 0, t. Therefore: 

j k j k 



la I.. ....I gt 

i7'i 
■Of.- 



(4.2.4) 



Thus in general to take the second derivative of a vertex one either adds two ^^s or an edge 



connecting to a trivalent vertex with two ^-s. The lowest-valence example is: 
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I 96 



I. .1?' 

7'j -''l': + 



(4.2.5) 



The second summand vanishes b is a based loop in each variable. As we will need it 

later, we record one further derivative: 



^".i, 9c 96 



I 1 9c I ' I I I I 

?0-,V. 9a 9c .1.96 9a ; I i 96 

\::v^::...-^ "-v-- ■■•■>^:- 



9a.i^ ..I.'?'' .;96 qa\^ ..1.'?'= ;96 




0^ 

Qo. . ..i.r I 96 



0+ +0 (4.2.6) 



4.3 Derivatives of the Green's function and the determinant 

In this section, we evaluate the derivatives of the Green's function G-y given by equation (2.1.10) 
with respect to qa-,t. We record also the derivatives of log 



dG 

It is easy to evaluate -a^- Let ^ be a based loop, and consider the derivative of equation (4.2.3): 



det 



dqodqi 



dq, 



+ 



(4.3.1) 



We contract with G^, which is a based loop in each variable, and use equations (4.2.1) and (4.2.3): 



.1 9a 



.9a 



.1. 9a 



(4.3.2) 



The first equality in equation (4.3.2) requires that is a based loop, which follows from differ- 
entiating G'y(?,0) = = G^{<^,t) with respect to qa- 

Evaluating ^ is not so easy, because it is not a based loop: from equation (4.3.9) it will follow 
that ^^^{t,t) = 6''^. Instead, we will differentiate equation (4.1.4) directly. To begin, we evaluate 
Q^Qq^ by differentiating equation (4.1.3) with respect to qa, using equations (4.2.4) and (4.2.5): 



I. 9a 

■91 q{--. 



,91 



:.7 



:.9l 



91 



9a I 



91 . 



+ ii) ^4-? + (if (B) + :;7 



9a. i 91 , ' ' 



91 



9a 



+ 



I. 91 



■7 :■ 



B) + 



. ■ / 



- V 91 



j,qa 

(1^1) + (7'; 



91 . ' " ~ . I 9a 



13 



9°. I., 91 



91 ^ ' ~ - ^ I 9a 



(4.3.3) 



We remark that \\B 



f ^ dB{qx) 
■•' ~ dqo 



We turn now to O^^I"'^'])"^- Recall equation (4.1.2). Then: 



1 I 

90 .V- ■■■■•■:!; 91 
X ■ 

t"'" 



90 .k-. ■■■■■■ ■■■■■■•■■■.1.91 90 .(■■•■■■'■ ■ ■ ' 

,■•■■ 91 ^' ^^91 ■■■•.. .•■•■■ 91 ^' 

■■; .1 .-. ■.'. ■■■■ .1 

-m) im) s + 



90...'^.-^^-^"^', ^ -.■■■■■•■•■••■',.91 



■■■•■■•■'■..91 



90 I 



10.,. ^.V"^\^.'. .1.91 



90 



91 



I 91 



(4.3.4) 



In the second line above 'we used the product rule and dropped any terms involving derivatives ■with 
respect to qo of functions only of qi. In the third line we used equations (4.1.3), (4.2.5) and (4.2.6). 
By the quotient rule: 



■■••■2.Tr~ rrT.i-- 



90.^-^^-':-'^9l 

K ■ 5 



90 



' 1 t 2[ 



•-.91 



91 



: 70 ■•. ..■ C7 ) 



(4.3.5) 



Thus, we can differentiate the Green's function by recalling equation (4.1.4) and using equa- 
tions (4.3.3) and (4.3.5). Many terms cancel: 



e(r-o 



(92[-5])-i 

1 -^0 

..J.. 91 90 J + 

ir 'r 



t 



\ 



•■I-TT" TTT-O- 

J 91 90 J 

;'7': (T- 



17~~ -^0 

90. -J.. 

^7'i 



+ J. 91 

(1: 



+ (? ^ r) 
(4.3.6) 
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Finally, consider the integrals implicit in the trivalent vertices in equation (4.3.8). Write p for the 
variable of integration for each, and consider the sum as a single integral /^_q- When p < ? < r, 
the second and third terms exactly cancel; when ? < p < t, the three terms agree up to sign; and 
when ? < r < /O, the first two terms cancel. (When t the contribution to the total expression is 
from the part we have abbreviate "(? -f-T- r)", and the analysis is similar.) But consider breaking the 

integral rsJ>^ as a similar sum of contributions from /> < ? < r, etc. Equation (4.1.4) implies 
that the contributions from each part exactly match the integrals in equation (4.3.8). This proves: 

^ 91 , - , 91 

q T q T <; r 

Given equations (4.2.6) and (4.3.4) and the well-known formula for the logarithmic derivative 
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of a determinant, we also have: 




(4.3.10) 
(4.3.11) 



We conclude this section by simplifying the right-hand side of equation (4.3.10) and the first term 
on the right-hand side of equation (4.3.11). We first expand the a-vertex: 




But the integrals implicit in the first two summands above are obviously symmetric under switching 
the two "7" legs. Recalling equation (4.1.4): 



Similarly: 




All together, we have the following simplification: 




(4.3.12) 



4.4 The cancellations 

In the previous sections, we evaluated the derivatives of every component of V^. In this section, we 
check equation (4.0.1), which in the graphical notation reads: 

■■S) = 2^^J --^J + 2 + --^J ® + 2 + 2 ® ® + © (^-^-l) 

We have adopted the following convention: the dashed lines (although not the dotted circles) count 
as edges and the new nodes count as vertices for the purposes of computing Euler characteristic, 
etc., and in equations with connected diagrams with different numbers of loops, we multiply each 
diagram T by {ih)^^^\ 

Recalling that = m + (y) log|det ^lagi^ I + X^{connected diagrams}/{symmetry}, the cal- 
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culations from the previous sections give: 



-& = '-^ + 2 ^:7>v^:^.^ + 2 I® } + 5: Autr ^'•'•'^ 

■ ■■■■■■■ r connected 

with one or 



In the sum, it is impHed that all other vertices are usual a-vertices and are of valence n > 3. The 

sum includes the diagram '■ 'j.i^ coming from ^ log|det ^g^gg^^ | • Otherwise, each diagram with a 

*.;7';- comes from the diagram formed by removing it (if the *.:7';- attaches to a trivalent vertex, its 
removal leaves just a solid edge), and the combinatorics are as in equation (4.1.1): if there are n 

equivalent ways to attach a *.:7';- to a diagram F, then |Aut ^ — r7,:| = ^ |Autr|. The second sum 

occurs similarly, with each '^V'y';^ _ .;7' replacing a Green's function, as per the second term in 
equation (4.3.9). 

Similarly, we have: 

^ (-^ff'^f) (4^4.3) 

■•..?:■■ ■■■■■ ^ Autr ^ ^ 

r connected 

11 

with one :'-y: 

(71 v . . .■■ (71 ^ ^ 

Then 2 — 2''.l''- ■■'7'- + a sum of diagrams each with one .;7';^ _ .;7";. , such that its removal 

yields a diagram that is still connected. 

The diagrams in ^ Vr,/] V-y f ^ come in two types. Other than the simply-connected diagram 

^ '^Vh''. im'^^ ! there are diagrams where the dashed line connects a complicated part to a ''V'b''.) 

and diagrams where the dashed line connects two non-simply-connected diagrams, which are not 
connected otherwise. But recall equation (4.1.3); then: 

1 1 9i>;.--'-'.:-*Ji 



2 ■& + ■& ^ + 2 



2* ^:-) + * ® + 2^7>v^7;+ E Autr ^'-'-'^ 

■ r connected 



with one /'-y'. or 
one ;"'y*. .•7*. 



But comparing equations (4.4.1), (4.4.2) and (4.4.4), we see that the only remaining terms are equal 
by equation (4.1.2). Thus we have proven Theorem 2.2.6. 

5 Proof of Theorem 2.3.3 

Recall that the pointwise convergence of distributions is defined by testing against compactly- 
supported functions. For a chosen open neighborhood O C M'^"'+^, let A be the set of all smooth 
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families of classical paths with boundary values in O, and fix {t,qi) G M""'"^. Throughout this 
section, we will let g : M"" — )■ M be smooth with compact support K such that {t} x K x {qi} C O. 
We will prove Theorem 2.3.3 in three steps. In Section 5.1, we characterize those families -f £ A 
such that J^''™^' exp(— (ifi)~^V^(t, qo, gi)) g{qo) dqo is non-zero, and show that there are only finitely 
many such families. In Section 5.2 we pick Oq C M" open with compact closure and show that 
for e > sufficiently small depending on Oq, for O = (0, e) x Oq x Oq there is a unique family 
of classical paths 7 with J^™^' exp(— (i^)~-^Ky(t, go, g{Qo)dqo 7^ 0. In Section 5.3 we compute 
the limit as t — >• of this formal integral. 

5.1 On the support of formal integrals 

We wish to understand the formal integral: 

formal 

exp{{ihy^Vy{t,qo,qi)) g{qo) dqo (5.1.1) 

s: 

We begin by splitting the exponent V-y = (—5'^) + {Vy + S^). Then V-y + S-y is valued in formal 
power series that begin in degree h, so exp((z/i)~^(Ky + 5.y)) is a formal power series. Thus the 
formal integral in equation (5.1.1) is determined by the ^ — )• asymptotics of integrals of the form 
/ K exp(-(z/i)-i5^(t, qo,qi)) g{qo) dqo. 

Such integrals are conventionally studied by the method of stationary phase. Indeed, by [4, 
Lemma 3.13], if 7^ in i^, then: 

/ exp{-{ih)-^Sy{t,qo,qi))giqo)dqo = Oih'^) ash^O (5.1.2) 
Jk 

dS 

By equation (2.1.8), ■Q^{t,qo,qi) = if and only if 7(t,go,gi)(0) = -^(^o), where 7(t,go,5i) the 
member of the family of paths 7 that connects go to qi in time t. If the family 7 contains a 
path with initial values (7(0), 7(0)) = {—B{qo),qo) for some qo G K, then it does so in an open 
neighborhood of qo, and 7 is determined on the connected component of {t} x K x {qi} containing 

d S 

{t,qo,qi)- Thus, as K is compact, there are only finitely many families 7 with -Q^{t,qo,qi) = 
for some qo G K, and so there are only finitely many paths 7 such that the formal integral in 
equation (5.1.1) is non-zero as a formal expression in h. 

5.2 Short-duration classical paths can be almost geodesic 

For e G M, consider the non-degenerate second-order ODE given by: 

Let (fe be the "flow for duration 1" map, i.e. ipe{v,q) = (7^(1), 7^(1)) for 7^ the unique solution 
to equation (5.2.1) with (7^(0), 7^(0)) = {v,q). It is defined on an open neighborhood of the 
section in TM" = M'^", it is smooth when it is defined, and it depends smoothly on e. Let </)e(f , q) = 
(g,7r o ip^[v,q)^, where vr : TM" — )• M" is the natural projection. Then (po is the isomorphism 
(v, q) I— >• {q,q + v). So the partial function (/)[_] : TM x M — )• M^" is defined on an open neighborhood 
of TR" X {0}. 

Fix Oo ^ with compact closure Oq, and let Oi D Oo be open with compact closure Oi. 
Let V = 00^(07 X Wi). It is a compact subset of TM" containing {0} x Oi. Then we can 
find eo > such that V x (— eo,eo) is contained in the domain of 4>[-]- Let {v,q) G V such 
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that q G Oi. For e S (— eo,eo), let 7^ be the duration-1 solution to equation (5.2.1) with initial 
conditions (7^(0), 7e(0)) = {v,q). Then 70 is nondegenerate, and the nondegeneracy condition 
depends smoothly on e. Thus for each (v, q) there is some number < ei{v, q) < eo so that for e G 
(— ei,ei), 7e is nondegenerate. By lemma 2.1.4, ei can be taken to depend lower-semicontinuously 
on {v,q). Thus it has a minimum value €2 on the compact set {{v,q) G V s.t. q G Oq}. 

Then €2 > satisfies the following: for each e G (—£2, £2) and for each {qo,qi) G Oq x Oq, we 
have chosen a nondegenerate duration-1 solution to equation (5.2.1) connecting go to qi, and the 
set of all these chosen paths for a given e is precisely the set of solutions 7^ to equation (5.2.1) with 
(7e(0),7e(0)) G V and 7(0), 7(1) G Oq. Thus this family 7[_] is smooth. 

However, if e > and 7^ is a duration-1 solution to equation (5.2.1), then r 1-^ 7e(e~^r) is 
a duration-e solution to equation (2.1.2) (equation (5.2.1) with e = 1). In particular, we have 
constructed a smooth family 7 of classical paths with boundary values ranging in Oq x Oq x (0, €2). 

Moreover, suppose that 7 is a duration-e solution to equation (2.1.2) with initial conditions 
7(0) G Oq and 7(0) = — i?(7(0)). Then if e is sufficiently small depending on 7(0), we have 
(e7(0),7(0)) G V, and so 7 is a member of our family. In particular, our family contains all paths 
with 7(0) = — i?(7(0)) that start and end in Oq with sufficiently small duration depending on the 
endpoints. 



5.3 The limit as t -> 

In Section 5.1, we showed that the only contributions to an integral of the form 



formal 
K 



exp((i?i) ^Vy{t,qo,qi)) g{qo)dqo 



formal 



K 



det 



92 [-5^ 



dqodqi 



exp(X;diagrams) g{qo) dqo (5.3.1) 



come from families 7 that include paths with initial conditions 7(0) = —-6(7(0)). Indeed, for fixed 
t,qi, equation (5.3.1) is supported only at the points qo £ K so that the duration-f path with 
initial conditions (7(0), 7(0)) = (— -B(go))Qo) ends at j{t) = qi. In Section 5.2 we showed that for 
O = (0, £2) X Co X Oq there is a unique such family, and that for qi G Oq and sufficiently small t 
depending on qi, there is a unique such point qQ. Indeed, we have qi — qo = — t-B(go) + O(t^). We 
remark also that exp(^ diagrams) = 1 + 0{K). 

We claim that for these triples {t,qo,qi), we have: 

1 \ kl 



(5.3.2) 



dqidqQ 

Indeed, consider integrating equation (2.1.2) twice, with initial conditions (7(0), 7(0)) = {vQ,qQ): 



^\t)=qi, + tvi,+ 
Then holding qQ fixed we have: 



^0 



dq^ dq^ 



dq^ 



d<, dt 



9=7(';) 



dvQ 



t5] 



JO 



dq^ \ dq 



+ 



dBi 
dq^ 

dBk 



dq^ dq'^ 




dc^dt (5.3.3) 



iJ=7(t) 
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We now take t to be sufficiently small and vq to vary only in a compact neighborhood of —B{qQ). 
Then ^ = 0(1) as t ^ 0, but then equation (5.3.3) gives = t^j + 0{t^). On the 

other hand, as we observed in equation (3.2.2), it follows from equation (2.1.8) that ( dqidqo ) ~ 

(It)"' = fe- This gives equation (5.3.2). 

We make a similar argument, now holding qi fixed: 



i ^ ,i /* r ffdBk dBA dC 



JO 



d-idt (5.3.4) 

q=7(?) 



^ = -t6i + 0{e) (5.3.5) 



„ 1 \ kli / -I \ Pel 



(ago)O ) = + ) ={i'''^ + 0{l))-'Y = t5^^ + 0{e) (5.3.6) 

Then, by [4, Theorem 3.15(ii)], we have: 

formal 

exp((i/i)"V^(t,go,gi)) 5'(^o)c?go = 
= (2^n)"/2 X e*W4 X det(t5i + 0(t2))-^/2 ^ ^-i^n)-^S-,{t,q,+B{q,)+0{t^),q,) ^ 

xdet{t6i+0{t^)Y^^ X g{qi + B{qi)+0{t^)) X (l + 0(?i)) 

But S,{t,q, + B{q,) + 0{t%q,) = + B{j) ■ 7 - ^(7)) dr = /o*(i(7 + B)^ + 0(1)) dr = 

/o (C>(t) + 0(1)) = 0(t). Thus: 

/formal 
exp{{ih)-X{t,qo,qi)) g{qo)dqo = {27: HT ' ^'"^^ ' \g {q^) + 0{h)) (5.3.7) 

What about the oidex-h terms? They can be given explicitly as a sum along the lines of our 
equation (2.2.7) — c.f. [14, 15, 8] — or implicitly as in [4]. All we need is the following fact: 
each degree in h is given by a finite sum in which for each summand there is some M so that the 
summand scales as { ^{Qq^yl^ ) times a product of strictly fewer than M terms given by derivatives 
of g and at qo = qi + 0{t^). 

The calculations in Section 4 in fact allow us to write down all derivatives of with respect 
to qo- In addition to the derivatives computed in Section 4.4, we have: 

{dqof- ^ _ |Aut(r,M)| ' -"^ ^^-"^-^^ 

connected marked 
Feynman diagrams (F, M) 

with n ordered .•'7'.s 

Consider expanding each edge as in equation (4.1.4). Then in each diagram a marked half edge 
converts some into = —B{qQ) + 0{t). On the other hand, each vertex corresponds to 

an integration J^* of bounded functions, and so the vertices are 0{t). Thus, by equation (5.3.2), 
J-^(r,M) = 0(tl^l+l^l) for each dia gram T in equation (5.3.8). As the number of diagrams at each 
order is finite, this proves that for A; > 3, = 0{t). 



20 



But the derivatives of g do not depend on t; = —-^ + ^{diagrams} = + ^{diagrams} = 
0{t); and (Q^y2 = "■^§^~p" + ^{diagrams} = 0{t^^). As each term after the first one in the asymp- 
totic expansion of equation (5.3.7) has strictly more occurrences of ( ^(igq')'^^ ) ^ than of derivatives 
of g and V^, we see that the 0(h) part of equation (5.3.7) is also 0{t). Thus, as formal expressions 
in h: 

/formal 
exp{{ih)-%{t, qo, qi)) g{qo) dqo = {2TThT'\'^^'^g[qi) (5.3.9) 

Theorem 2.3.3 is proved. 
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